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Methanogens typically occur in reduced anoxic environments. However, in recent studies 
it has been shown that many aerated upland soils, including desert soils also host 
active methanogens. Here we show that soil samples from high-altitude cold deserts 
in the western Himalayas (Ladakh, India) produce CH4 after incubation as slurry under 
anoxic conditions at rates comparable to those of hot desert soils. Samples of matured 
soil from three different vegetation belts (arid, steppe, and subnival) were compared 
with younger soils originating from frontal and lateral moraines of receding glaciers. 
While methanogenic rates were higher in the samples from matured soils, CH4 was 
also produced in the samples from the recently deglaciated moraines. In both young 
and matured soils, those covered by a biological soil crust (biocrust) were more active 
than their bare counterparts. Isotopic analysis showed that in both cases CH4 was 
initially produced from H2/CO2 but later mostly from acetate. Analysis of the archaeal 
community in the in situ soil samples revealed a clear dominance of sequences related to 
Thaumarchaeota, while the methanogenic community comprised only a minor fraction of 
the archaeal community. Similar to other aerated soils, the methanogenic community was 
comprised almost solely of the genera Methanosarcina and Methianocella, and possibly 
also Methanobacterium in some cases. Nevertheless, ~10'^ gdw^^ soil methanogens 
were already present in the young moraine soil together with cyanobacteria. Our results 
demonstrate that Methanosarcina and Metlianocella not only tolerate atmospheric oxygen 
but are also able to survive in these harsh cold environments. Their occurrence in 
newly deglaciated soils shows that they are early colonizers of desert soils, similar to 
cyanobacteria, and may play a role in the development of desert biocrusts. 
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INTRODUCTION 

Methanogenic archaea are strict anaerobes and are therefore 
typically found in permanently water-logged, anoxic habi- 
tats such as wetlands, rice fields, marine sediments, hot 
springs, and the guts of ruminants and termites (Zinder, 1993; 
Liu and Whitman, 2008; Brune, 2011). Despite that, it has 
been shown that active methanogens inhabit many aerated, 
oxic soils including desert soils (Peters and Conrad, 1995; 
West and Schmidt, 2002; Angel et al., 2012). Methanogens 
belong to the phylum Euryarchaeota and are phylogeneticaUy 
divided into six orders: Methanobacteriales, Methanococcales, 
Methanomicrobiales, Methanosarcinales, Methanopyrales (Ferry, 
1994), and Methanocellales (Sakai et al, 2008). Recently two 
archaeal strains from the order Thermoplasmatales isolated from 
human and termite guts have been shown to be methanogenic 
and were reclassified as the seventh order of methanogens — the 
Methanoplasmatales (Dridi et al, 2012; Paul et al, 2012). WhUe 



the aforementioned anoxic environments typically host a wide 
variety of methanogens, only the genera Methanosarcina and 
Methanocella were found in aerated soils, and it has been sug- 
gested that these methanogens are universal inhabitants of upland 
soils (Angel et al., 2012). All known methanogens possess the gene 
encoding for the enzyme methyl coenzyme M reductase, which is 
responsible for the last step of methanogenesis in all known path- 
ways. Thanks to its universal occurrence and conserved sequence, 
the gene encoding for its a-subunit — the mcrA — is commonly 
used as a phylogenetic marker gene for methanogens (Lueders 
et al, 2001; Friedrich, 2005). 

The production of biogenic CH4 is important in nature since 
it is the terminal step in the biodegradation of organic matter 
under anoxic conditions (Deppenmeier et al., 1996). In most ter- 
restrial environments CH4 arises primarily from the reduction 
of CO2 (hydrogenotrophic methanogenesis) or the cleavage of 
acetate (aceticlastic methanogenesis). Although the two processes 



www.frontiersin.org 



December 2013 | Volume 4 | Article 359 | 1 



Aschenbach et al. 



Methanogens in high-altitude cold desert soils 



yield the same product, the source of the CH4 can nevertheless 
be differentiated by analysing the stable isotopic signatures of the 
carbon in the CH4 and its precursors (acetate and CO2; Conrad, 
2005). The production of biogenic CH4 in nature typically leads 
to its release to the atmosphere where it acts as a greenhouse 
gas which is 25-33 times more potent than CO2 (Shindell et al., 
2009). 

Deserts (non-polar arid and hyperarid regions) are the 
largest biome on Earth, spanning over 20% of the land sur- 
face (Middleton and Thomas, 1997). Many of these deserts are 
characterized by a warm climate but some, particularly in high 
altitudes, experience frequent sub-zero temperatures through- 
out the year, and are termed cold deserts. The flora, fauna, and 
microbiota in these regions are therefore adapted to coping with 
low temperatures in addition to the desiccation stress common 
to all deserts (Dvorsky et al., 2013). The area of East Ladakh, 
India in the western Himalayan slopes is an extensive and thinly 
populated high-altitude plateau characterized by both low tem- 
peratures as well as low precipitation owing to its location in the 
Himalayan rain shadow (Dvorsky et al, 2011). This area can thus 
serve as a nearly pristine model for studying the ecology of cold 
deserts. 

The soils of arid and semi-arid regions are often covered by 
a unique layer, a few millimeters thick and densely colonized by 
microorganisms, termed biological soil crust (biocrust). Biocrusts 
are formed by living organisms and their by-products, creating a 
layer of soil particles bound together by organic materials, e.g., 
the sticky sheath material of Microcoleus (Biidel, 2003). Biocrusts 
are predominantly composed of bacteria (cyanobacteria and oth- 
ers), archaea, fungi, and green and brown algae (Belnap et al., 
2003; Angel and Conrad, 2013). In more humid areas, mosses, 
lichens, and liverworts can also be present. The particular com- 
position of microorganisms that dominate the crust varies in 
nature and is largely determined by climate, soil type, the suc- 
cessional stage of the crust development as well as by historical 
contingency (Biidel, 2003; Bahl et al., 2011). In contrast to hot 
desert biocrusts, which are smooth and very thin, in cold climates 
biocrusts tend to be much thicker and develop three-dimensional 
structures, up to several centimeters high, in response to freeze- 
thaw cycles, termed as rolling and pinnacled crusts (Belnap, 
2003). Whether this morphological difference is also reflected in 
differences in microbial community composition and function 
still requires research. In previous studies it has been shown that 
biocrusts from hot deserts can produce CH4 when wet (Angel 
et al, 2012). The primary goal of this study was therefore to 
investigate whether CH4 production can also be observed in high- 
altitude cold deserts such as in Eastern Ladakh. In addition, the 
presence of receding glaciers allowed us to also study how quickly 
methanogens colonize young soils. We hypothesized that similar 
to marine microbial mats (Hoehler et al, 2001), methanogens 
are an integral part of biocrusts and will be found already at 
early successional stages of the soil development together with 
cyanobacteria. Since it has been estimated that ~40% of the land 
area in Ladakh is covered by biocrusts (Janatkova et al., 2013), 
we also set out to compare the potential activity and commu- 
nity composition between patches covered by a biocrust and 
bare soil. 



For this purpose, soil and biocrust samples from three different 
vegetation belts: arid, steppe, and subnival, as well as from front 
and lateral moraines of receding glaciers were collected and tested 
for the production of CH4. In addition, the carbon isotopic sig- 
natures of CH4 and CO2 were analyzed to determine the pathway 
by which CH4 was formed. Finally, the methanogenic community 
was characterized by molecular analysis. 

MATERIALS AND METHODS 

SOIL SAMPLES AND SITE CHARACTERIZATION 

Soil samples were collected from two sites in a high-altitude cold 
desert in Ladakh, India, from Nubra Valley and from Tsomoriri 
Plains. Nubra Valley is characterized by a dry climate (100 mm 
annual precipitation) and a neutral to alkaline soil pH. In Nubra 
Valley samples were collected from the frontal moraines of three 
glaciers (Nl - N34° 40' 26.87" E77° 45' 36.17"; N2 - N34° 
39' 16.98" E77° 45' 8.63" and N3 - N34° 39' 16.10" E77° 
44' 36.48") located at 5400, 5300, and 5150 m above sea level, 
respectively. Three to six frontal moraines where studied at each 
glacier, and from each moraine three samples were collected at 
different positions (A, B, and C; Figure Al). The annual mean 
temperatures in these sites range from —1.6 to — 3.6°C (mea- 
sured at 4850 and 5250 m between Aug. 2009 and Aug. 2011; 
Dvorsky et al, 2013). The second sampling site was Tsomoriri 
Plains where the annual precipitation is just under 100 mm 
yr^' and the soil pH is neutral to alkaline. At this site, soil 
and crust samples from three vegetation belts were sampled: 
arid, steppe, and subnival located at 4700, 5300, and 5800 m, 
respectively (N 32°58'56,51" E78°2l'24,95", N 32°59'31,75", 
E 78°24'7,56", N 33° O' 23,26" E 78° 26' 46,76"). Furthermore, 
samples from lateral moraines of the Chamser glacier (N32° 
59' 13.17" E78° 25' 54.084") were also collected at altitudes 
5650, 5700, and 5800 m. Here, samples were collected from lat- 
eral moraines (A, B, and C), each divided into three transects 
(TSl, TS2, and TS3; Figure Al). From each moraine/transect 
samples were collected from the top and the bottom of the 
moraine. The annual mean temperatures at this site range from 
-4.4 to -10.4°C (measured at 5350 and 5850 m between Aug. 
2009 and Aug. 2011; Dvorsky et al, 2013). Sample character- 
istics are given in Table 1. The sampling areas were partially 
covered by a biocrust, and soil samples were collected from 
either the top soil which was covered by biocrust or from bare 
areas (0-5 cm), which were in proximity. In Nubra Valley only 
biocrust samples were collected. Samples of approximately 100- 
150 g soil from the top layer were collected into Whirl- Pak® 
bags (118-ml sterile sampling bags, Nasco) and air dried in 
the field. Altogether a total of 84 samples were collected, and 
shipped at room temperature to Germany for further analy- 
sis. For nucleic acids extraction, 1 g of each sample was split 
in two 2 ml tubes upon arrival at the laboratory and stored 
at — 80°C until extraction. Soil pH was determined in a 1:1 
soihwater slurry solution. Soil content of total and organic car- 
bon and total nitrogen were analyzed using an elemental analyser 
(vario Micro cube, CHNS mode, Elementar) by the Analytical 
Chemical Laboratory of the Philipps-Universitat, Marburg while 
the stable carbon isotope signature (S'^C) of the total and 
organic carbon were analyzed at the Institute for Soil Science 
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Table 1 | Physico-chemical characteristics and methanogenic potentials of a selection of the soil samples. 



Location 


Site 


Analyzed sample 


pH 




Ctot (%)^ 


Corg (%)' 


Carbonate (%) 


Lag [d] 


CH4 rate'' 




















(nmol g~^ d~M 


Nubra 


Glacier 1 


N1/3A Crust 


8.6 


0.12 ±0.01 


2.82 ± 0.02 


1.19 ±0.03 


58.0 


31 


169.8 ±81 .7 


Valley 


Glacier 2 


N2/4A Crust 


8.7 


0.07 ±0.01 


1.08 ±0.08 


1.02 ±0.24 


4.9 


22 


4.8 ±0.5 




Glacier 3 


N3/1C Crust 


8.2 


0.12 ±0.01 


2.46 ±0.10 


0.88 ±0.03 


64.2 


13 


68.4 ±43.0 


Tsomoriri 


Chamser 


TS 


7.5 


0.24 ±0.01 


2.5 ± 0.14 


2.53 ± 0.17 


BLD 


18 


849.8 ± 143.2 


Plains 


Glacier 


3B-bottom/crust 




















TS 3B-top/bare soil 


7.5 


BLD'^ 


0.59 ±0.03 


0.56 ±0.02 


4.3 




BLD 




Arid 


Crust 1 


8.2 


0.03 ±0.01 


0.63 ±0.07 


0.28 ±0.01 


55.8 


16 


378.0 ± 185.6 






Bare soil 1 


8.7 


0.01 ±0.00 


0.16 ±0.05 


0.14 ±0.01 


13.8 


21 


8.9 ±2.4 




Steppe 


Crust 2 


7.7 


0.10 ±0.02 


0.81 ±0.06 


0.81 ±0.00 


13.5 


16 


21 .8 ±7.7 






Bare soil 2 


8.0 


0.02 ±0.00 


0.38 ±0.04 


0.25 ±0.00 


34.0 


19 


19.9 ± 16.8 




Subnival 


Crust 3 


7.0 


BLD 


0.96 ±0.11 


0.83 ±0.08 


0.6 


16 


420.1 ±27.2 






Bare soil 3 


6.9 


0.02 ±0.00 


0.30 ±0.04 


0.29 ±0.05 


3.1 


21 


57.9± 13.9 



"Samples were analyzed In duplicates. 

^ Means of the 3 incubated technical replicates. 

"BLD, below the limit of detection. 



and Forest Nutrition (IBW) at the University of Gottingen, 
Gottingen, Germany, using an elemental analyser (NA2500, CE- 
Instruments) coupled over an interface (Conflo III, Thermo) to 
a mass spectrometer (Delta plus, Finnigan MAT). In both cases 
measurements were done before and after acidification with 10% 
HCl, and the difference in values was attributed to carbonate 
(Nusslein et al, 2003). 

INCUBATION CONDITIONS AND GAS MEASUREMENTS 

For determining the methanogenic potential of a sample, 5 g 
of sieved soil were incubated in a 27-ml pressure tube and 
amended with 5 ml sterile, distilled-deionized water in triplicates. 
The tubes were closed with butyl rubber stoppers (cleaned by 
boiling before use), purged with N2 and incubated at 25°C in 
the dark. Samples were sacrificed after CH4 concentration has 
reached about 15,000 nmol g- dry- weight^' (gdw) or after 95 days 
if less CH4 was produced. Previous experience suggested that 
such a level of accumulated CH4 should allow for easy detection 
of methanogens through TRFLP (data not shown). Each sam- 
ple was incubated in triplicated tubes representing 3 technical 
replicates. Prior to gas analysis, the tubes were shortly shaken 
by hand to equilibrate gas and aqueous phase. Measurements of 
CO2 and CH4 were performed following Angel et al. (2012), gas 
samples (500 |il) were taken from the headspace using a 500-(jl1 
glass gas-tight pressure-lock syringe (Vici) every 7-14 days and 
analyzed immediately. Methane and CO2 concentrations were 
analyzed using a GC (GC-8A; Shimadzu) equipped a 3 m, 0 1/8" 
stainless steel column filled with Hayesep Q 80/100 mesh, a meth- 
anizer (Ni-catalyst at 350°C, Chrompack) and a flame ionization 
detector (SRI, temperature: 160°C). The injector temperature 
was 160°C, the oven temperature: 120°C, and the carrier gas 
was H2. 

The isotopic signatures of the carbon in the CH4 and CO2 were 
determined using a Gas Chromatograph Combustion Isotope 
Ratio Mass Spectrometer (GC-C-IRMS; Trace GC Ultra, Thermo 



Fischer Scientific), following Conrad et al. (2009). The prin- 
ciple operation of GC-C-IRMS has been described by Brand 
(1996). The CH4 and CO2 in a 500 |il gas samples were first 
separated in a Trace GC Ultra Gas Chromatograph using a 
Pora PLOT Q column (27.5 m length, 0.32 mm i.d.; 10 [im film 
thickness; Varian, Palo Alto, CA, USA) at 30° C with helium 
(99.996% purity; 2.6ml/min) as carrier gas. After conversion 
of CH4 to CO2 in a GC Isolink 1030 Oxidation Reactor 
at 940°C, the isotope ratio of ^^C/'^C was analyzed in an 
IRMS (Delta V Advantage, Thermo). The isotope reference gas 
was CO2 (99.998% purity; Air Liquide), calibrated with the 
working standard methylstearate (Merck). Values are reported 
in the delta notation against the Vienna Pee Dee Belemnite 
(Hayes, 1993): 

S'^C = W^RsJRst - 1) (1) 

Where R^^ ="CJ^^Cs^ and Rst =^^Cst/^^Cst of sample (sa) and 
standard (st), respectively. 

The preference of a reaction for the light carbon isotope is 
described by the fractionation factor (a), which is defined as: 

a = (Ssubstrat + 1000)/(8product + 1000) (2) 

equivalent to: 

E = (l-a)1000 (3) 

Where Ssubstrat and Sproduct are the S'^C values of the substrate 
and the product, respectively. 

The relative fraction of the hydrogenotrophicaUy derived CH4 
was calculated with the following equation: 

/h2 = (S^^CcH4 — Sma)/(Smc — Sma) (4) 
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Where 8'^Cch4 = S'^C of CH4 in the headspace; is the S^^C 
of CH4 produced solely from acetate and Smc is the S^^C of CH4 
produced solely from H2/CO2. Values of 8ma and Smc were esti- 
mated from measured 8^^Cco2 (S^^C of CO2 in the headspace) 
using fractionation factors (a or e values) from the literature and 
assuming that the 8^^C of acetate (S^^Cac) is identical to the S^^C 
of organic matter (S^^Corg): 

Sma = S^^Corg + Sacetate,CH4 (5) 
8mc = SC02 + 6C02,CH4 (6) 

EXTRACTION AND AMPLIFICATION OF NUCLEIC ACIDS 

After incubations were completed, the pressure tubes were 
opened and the slurries were centrifuged at 10,000 rpm for 2 min 
to remove pore water. One gram of soil slurry was frozen in liquid 
nitrogen and stored at — 80°C. Total nucleic acids were extracted 
from 0.5 g of the dry unincubated soil as well as from the incu- 
bated slurry samples as previously described (Angel, 2012; Angel 
etal, 2012). 

The primers for the PGR reactions are listed in Table 2. 
Each PGR reaction was 50 \l\ in volume and contained: 10 (xl 
GoTaq®Flexi 5x Green Buffer (Promega), 0.2 mM dNTP mix- 
ture, 1.5 mM MgGl2, 0.8|xg/ixl BSA (Roche), 0.25 (jlM of each 
primer, 1.5 U of GoTaq® DNA polymerase (Promega), and 1 |xl 
of DNA template. The following programme was used: 94° G for 
4 min followed by 30 cycles of 94° G for 30 s, 52° G for 30 s, and 
72° G for 45 s and a single step of final elongation at 72° G for 
10 min. 

For analysis of terminal restriction fragment length poly- 
morphism (TRFLP) of archaeal 16S rRNA genes, the forward 
primer 109f was labeled with 6-FAM (6-carboxyfluorescein) at 
the 5' end. Approximately 200 ng of purified PGR products 
were digested overnight at 65° G with the restriction enzyme 
Taq"l (New England BioLabs). TRFLP sample preparation and 
processing were performed as previously described (Angel et al. 



2012). For analysis, the height of the measured peaks was used, 
and expressed in relative abundance of the total height. All peaks 
under 1% of the total height were treated as noise and were 
removed. 

CLONING AND SEQUENCING AND PHYLOGENETIC ANALYSIS 

Four samples (TS 2B-top/biocrust, TS 3B-bottom/biocrust, 
subnival-biocrust 3, and subnival-bare soil 2) showing the most 
diverse TRFLP patterns were used for cloning of archaeal 16S 
rRNA genes and Sanger sequencing, in order to assign TRFs 
to genus. Libraries were constructed using purified PGR prod- 
ucts (GenElute PGR cleanup kit, Sigma) which were cloned into 
the pGEM-T easy vector (Promega) following the manufacturer's 
instructions. Sanger sequencing services were provided by GATG 
(Germany) and M13f and M13r primers were used to sequence 
from both flanking regions of the vector to the insert. Twenty 
four clones from each sample were analyzed (96 clones in total). 
Only reads that fulfilled the following criteria were used: the two 
reads were overlapping up to the opposite primer sequence, the 
assembled contig was longer than 800 bp, the contig was con- 
firmed to be archaeal and also not chimeric using BLAST. Sixty six 
sequences passed our quality filtering and were used for phyloge- 
netic analysis. Sequences were aligned against the SILVA 108 SSU 
Ref database (Quast et al., 2013) using the ARB software pack- 
age (Ludwig et al, 2004). The phylogenetic tree was calculated 
using RAxML (Stamatakis, 2006), implemented in Arb, using 
rapid hill climbing algorithm, PROTMIXJTT evolutionary model 
and a 1000 bootstrap runs. Nucleotide sequences were submit- 
ted to GenBank and can be found under the following accession 
numbers: KF445438-KF445503. 

QUANTITATIVE REAL-TIME PCR ASSAYS 

Quantitative real-time PGR (qPGR) assays were used to quan- 
tify the 16S rRNA gene of cyanobacteria and the mcrA gene of 
methanogens. The assays were based on SYBR® Green and were 
performed as previously described (Angel et al., 2011; Angel and 



Table 2 | Overview of the used primers in this study. 



Oligo. 


Target 


Oligo. sequence {5'-3') 


Position'' 


GC (%) 


Tm 


Amplicon 


Assay 


References 


name'' 












size 






ARC 109-F 


Archaea 16S rRNA gene 


ACK GOT GAG TAA CAC GT 


109-125 


47 


54 


826 


TRFLP 
sequencing 


GrolSkopf 
et al., 1998 


ARC 934-R 




GTG GTG CCC CGC CAA TTC CT 


91 5-934 


65 


71 








CYA 


Cyanobacteria 16S rRNA 


GRG GAA TYT TCC GCA ATG GG 


359-378 


60 


63 


447 


qPGR 


Nuebel 


359_mod-F 


gene 














etal., 1997 


CYA 781-R 




GAG TAG WGG GGT ATG TAA TCC 
CWTT 


791-805 


46 


64 








mlas-mod-F 


Universal mcrA gene 


GGY GGT GTM GGD TTG ACM 
CARTA 


976-998a 


43-65 


68 


469 


qPCR 


Steinberg 
and Regan, 
2008 


mcrA-rev-R 




GGT TCA TBG GGT AGT TVG GRT 
AGT 


1421-1444a 


42-54 


66 









forward primer; R, reverse primer 

''Positions are based on tfie foilowing: primers targeting the IBS rRNA gene — E.coli; primers targeting the mcrA gene — M. thermautotrophicus mcrA gene accession 
number: U 10036 (following Steinberg and Regan, 2008). 
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Conrad, 2013). A standard containing a known number of DNA 
copies of the target gene was used for all assays in serial dilutions 
for generating a calibration curves. These standards were gener- 
ated from an environmental clone of a cyanobacterium and from 
a pure culture of Methanosarcina thermophila. All qPCR reac- 
tions were performed on an iCycler thermocycler equipped with 
MyiQ™ detection system (BioRad). The resulting data were ana- 
lyzed with the iQ5 optical systems software (BioRad). Accounting 
for the dilutions of DNA template, our assays allowed detecting 
5 X 10^ copies per gram dry soil. 

RESULTS 

SOIL CHARACTERISTICS 

The pH in the samples obtained from the vegetation belts and the 
lateral moraines (Tsomoriri) was neutral to alkaline (pH 6.93- 
8.66), which is typical for desert soils, while the samples from 
the frontal moraines of Nubra Valley were slightly more alkaline 
(pH 8.14-8.96). Overall the samples were poor in both carbon 
and nitrogen. Total carbon content ranged from 0.16 to 2.8% 
while that of nitrogen was below the detection limit in many sam- 
ples and reached at most 0.24% of the dry soil mass (Table 1). 
Further analysis showed that in many of the samples much of 
the carbon was inorganic, i.e., carbonate, which made up as 
much as ~60% of the total carbon in the samples (Table 1). 
The isotopic signature of the carbon in the soil typically ranged 
from —10 to — 22%c and tended to be heavier in samples with 
higher carbonate content (Table 3). 

METHANE PRODUCTION POTENTIAL AND THE STABLE ISOTOPE RATIO 
OFCH4 ANDCO2 

Methane production was tested in all 84 samples. The CH4 pro- 
duction rate was always higher in the biocrust samples than in 
the bare soil samples. In total, 50 (60%) of the samples produced 
at least some CH4 ( > 1 nmol gdw^ ^ d^ ' ; Table 1 ; Table Al ) . The 
CH4 production rate of 14 highly methanogenic samples (~30% 
of the methanogenic samples) reached over 100 (and in some 



cases up to around 1000) nmol gdw^^ d^^ Linear CH4 accu- 
mulation rates were detected in the vegetation belt samples at 
Tsomoriri Plains after a lag phase of 19 days on average (Table 1; 
Table Al). First traces of CH4 were detected in the lateral moraine 
samples at Tsomoriri Plains after a lag phase of 24 days on average 
and in Nubra Valley samples after 3 1 days on average. 

All samples from the different vegetation belts produced CH4 
(Table 1 ; Table Al ) . The highest methanogenic activity was found 
in the subnival soil samples with 835 ± 228 and 65 ± 5 nmol 
CH4 gdw^^ d^^ on average in the crust and bare soil, respec- 
tively, followed by the arid samples with 329 ±124 and 19 ± 
14 nmol CH4 gdw^^ d^' on average in the crust and bare soU, 
respectively and by the steppe samples with only 99 ± 69 and 17 
± 6 nmol CH4 gdw^^ d^^ on average in the crust and bare soU, 
respectively. 

Also the soils from the lateral moraines of the receding 
Chamser glacier in Tsomoriri Plains showed methanogenic 
potential, though not in all samples and in lower amounts com- 
pared to the vegetation belts samples. In transect 1 samples did 
not produce CH4, in transect 2, only 4 of the 10 samples produced 
CH4, and in the last transect (3), 9 of the 10 samples produced 
CH4. Here as well, biocrust samples produced more CH4 than 
the bare soils in nearly all cases except moraine C, where the 
bare soU sample produced slightly more CH4 than the sample 
with the crust (Table Al). On average, the active biocrust sam- 
ples produced CH4 at a rate of 239 ± 114 nmol CH4 gdw^^ d^^ 
while the active bare soil samples produced only 89 ± 48 nmol 
CH4 gdw-i d-i. 

In Nubra Valley a larger proportion of the samples were active 
compared to the Tsomoriri Chamser glacier samples and all three 
glacier moraines showed at least some CH4 production, but most 
samples produced very little CH4. On average, active samples 
here (which were all biocrust) produced 20 ± 8 nmol CH4 gdw^' 
d-i. 

In addition to the measurement of gas concentrations, the iso- 
topic signatures of the carbon containing compounds, namely 



Table 3 | Stable carbon isotope analysis. 



Location 


Site 


Analyzed sample 


813Ctot 

(%o) 


S13Corg 

(%o) 


S13CcH4 

Begin''''' 


813CcH4 

End 


513Cco2 
Begin^'' 


813Cco2 
End 


fHzBegin^'' 


fH2End 


Nubra Valley 


Glacier 1 


N1/3A Crust 


-10.3 


-23.6 


-71 


-38.7 


-18.4 


-20.7 


0.7-1.0 


0.1-0.2 




Glacier 2 


N2/4A Crust 


-20.5 


-25.2 


-66.4 


-50.4 


-24.1 


-24.1 


0.6-0.9 


0.3-0.5 




Glacier 3 


N3/1C Crust 


-10.5 


-23.1 


-58.2 


-53.8 


-21.7 


-20.4 


0.4-0.7 


0.4-0.6 


Tsomoriri 


Chamser 


TS 3B-bottom/crust 


-22.2 


-22.4 


-64.1 


-39.8 


-14.6 


-19.1 


0.5-0.8 


0.1-0.2 


Plains 


Glacier 


TS 3B-top/bare soil 


-20.9 


-21.2 


BLD= 


BLD 


BLD 


BLD 


BLD 


BLD 




Arid 


Crust 1 


-14.0 


-19.7 


-52.7 


-47.4 


-11.9 


-14.6 


0.4-0.5 


0.3-0.4 






Bare soil 1 


-13.6 


-24.0 


-61 


-64.2 


-24.1 


-22 


0.5-0.8 


0.6-0.8 




Steppe 


Crust 2 


-17.5 


-20.4 


-69.7 


-70 


-16.5 


-19.7 


0.6-1.0 


0.6-1.0 






Bare soil2 


-22.3 


-23.9 


-69.6 


-52.4 


-23.6 


-22.8 


0.6-1.0 


0.3-0.5 




Subnival 


Crust 3 


-21.6 


-21.3 


-69.8 


-43.5 


-4.4 


-12.4 


0.6-1.0 


0.2-0.3 






Bare soil 3 


-22.3 


-22.1 


-62.6 


-41.4 


-27.0 


-22.4 


0.5-0.8 


ca. 0.2 



^ Means of the 3 technical replicates. 

^ "Begin " and "end" refer to the first and last measurements. 
'^BLD, below the limit of detection. 
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CH4 and CO2, were measured in this study. The isotopic signa- 
tures (8^-''C) of the carbon in the CH4 and CO2 indicate which 
methanogenic pathway — hydrogenotrophic or aceticlastic — was 
favored in our incubations. The 8^^C of the CH4 ranged 
from — 100 to — 30%c, and showed generally a positive trend with 
time (and therefore also with concentration; Figures 1-3). At the 
beginning of the incubations, the average 8^^Cch4 values of all 
analyzed samples were around —70 to — 52%c (Figures 1-3). Over 
time values shifted to heavier 8^^C in the soil samples from the 
lateral (8^^Cch4 = — 40%c; Figure 2) and the frontal moraines 
(8^^CcH4 = —54 to — 32%c; Figure 3). In contrast, the isotopic 
signature of CH4 of the steppe-biocrust sample as well as the 
bare soil samples from arid sites were relatively constant over time 
(8^^CcH4 = -70 to 64%c; Figure 1). As for CO2, the h^^C ranged 
from —30 to — 5%c (Figures 1-3). The bare soil samples showed a 
heavier isotopic signature in the CO2 compared with the biocrust 
samples. 

The characteristic S'^C values of CH4 and CO2 at the begin- 
ning and the end of incubation are summarized in Table 3. These 
values were used to calculate the fraction of hydrogenotrophic 
methanogenesis {/hi) using Equation (4). For these calcu- 
lations we assumed Eac,cH4 = — 10%c (Goevert and Conrad, 
2009; Angel et al, 2012) and ?>^^C^c = 8^^Corg- On average, 
Sma was therefore — 32.2%o. We further assumed previously 
published fractionation factors for conversion of H2/CO2 to 
CH4 (-73 and -49%o; Fey et al, 2004; Conrad et al., 2010) 
so that 8inc was in the range of —90 to — 70%c. The rela- 
tive fraction of CH4 derived from H2/CO2 in the early emit- 
ted methane was relatively high [fni > 0.50 in almost all 



samples; Table 3). Thereby the samples seem to be generally 
dominated by hydrogenotrophic methanogenesis in the begin- 
ning. Toward the end of the incubation, however, fyi dropped 
to 0.10-<0.50. Nevertheless, in some samples (arid-bare soil, 
steppe-biocrust, and glacier 3-lC biocrust) hydrogenotrophi- 
cally derived CH4 remained a dominant pathway over time 
(/h2 = 0.4-1.0). 

COMPOSITION OF THE ARCHAEAL COMMUNITY 

TRFLP was used to characterize the archaeal community in dif- 
ferent samples. TRFLP profiles were generated only for a subset 
of the total samples, i.e., for one original sample (before incuba- 
tion) from each vegetation belt and glacial moraine at Tsomoriri 
or Nubra, and for two out of three of the respective incu- 
bation technical replicates. The tested samples showed similar 
occurrence of TRFs but their proportion varied between sites 
(Figure 4). Altogether some 14 different TRFs could be reliably 
detected in the analyzed samples, but only five of the major 
TRFs could be identified in the sequences from the clone library 
(Figure A2). 

The diversity of archaea in the in situ samples was domi- 
nated (80-100% of TRFs) by Thaumarchaeota (TRFs 530 and 
602 bp) while no methanogens could be detected. Following incu- 
bation, the diversity of archaea increased with the appearance of 
methanogens. Methanosarcina (608 TRF) and Methanocella (295 
TRF) appeared in all CH4 producing samples and in one case also 
Methanobacterium (702 TRF) was detected. The proportion of 
methanogens in the incubated methanogenic samples increased 
to 30-100% of the total community, with the level of enrichment 
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FIGURE 1 I Values of s"C of CH4 and CO2 as a function of CH4 and CO2 
concentrations in incubations with soil samples from the different 
vegetation belts In Tsomoriri Plains. Depicted are the three vegetation 
belts samples: arid, steppe, and subnival. (A) S^'^C of CH4 and (B) S^'^C of 



CO2. Only samples producing >40 nmol g"^ CH4 produced reliable S^-'C 
measurements of CH4 and are shown here. Temporal progression is always 
along the x-axis. Three biological replicates were analyzed and each data point 
represents a mean of n = 3 technical replicates. 
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FIGURE 2 1 Values of of CH4 and CO2 as a function of 

CH4 and CO2 concentrations in incubations with soil samples 
from the different lateral moraines in Tsomoriri Plains. (A) S^'^C 
of CH4 and (B) S^-^C of COa- Only samples producing >40nmol 



g"' CH4 produced reliable 8"C measurements of CH4 and are 
shown here. Therefore, only transects TS2 and TS3 are plotted. 
Temporal progression is always along the x-axis. Means of n = 3 
(technical replicates). 
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FIGURE 3 I Values of 8^^C of CH4 and CO2 as a function of CH4 and 
CO2 concentrations in incubations with soil samples from the 
different frontal moraines in Nubra Valley. (A) S'^C of CH4 and (B) S'^C 
of CO2. Only samples producing >40nmol g^' CH4 produced reliable S'-'C 



measurements of CH4 and are shown here. Nl, glacier 1 in Nubra Valley, 
containing 4 moraines; N2, glacier 2 with 6 moraines, N3, glacier 3 with 4 
moraines. Temporal progression is always along the x-axis. Means of n = 3 
(technical replicates). 



correlating with the methanogenic potential. The nine TRFs (63, 
71, 99, 105, 395, 507, 574, 581, 661 bp) which could not be 
affiliated to the sequence data occurred only occasionally in the 
samples, and made up only ~14% of the population. In the case 



of one sample from Nubra Valley (N2/4A), unidentified TRF = 
99 bp was shown to contribute 27% to total T-RFs. However, the 
CH4 production rate in this sample was very low, thus this TRF is 
unlikely to have represented a methanogen. 
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FIGURE 4 I TRFLP community profiles of archaeal 16S rRNA genes 
from different soil samples. The community composition was 
characterized in a selection of in situ and incubated soil samples 
from the three vegetation belts, the lateral moraines in Tsomoriri 
Plains and the frontal moraines in Nubra Valley. For the incubated 



samples technical replicates were analyzed. Means ± SE of n = 2 
(technical replicates). TS1, Tsomoriri moraine 1; TS2, Tsomoriri 
moraine 2; TS3, Tsomoriri moraine 3; N1/3A, Nubra glacier 1 
moraine 3a; N2/4A, Nubra glacier 2 moraine 4a; N3/4B, Nubra glacier 
3 moraine 4b. 



QUANTIFICATION OF CYANOBACTERIA AND METHANOGENS ALONG A 
LATERAL MORAINE 

To determine how quickly methanogens colonize soils that were 
newly exposed by deglaciation and to better understand the rela- 
tionship between methanogens and biocrusts, we quantified the 
density of methanogens and cyanobacteria in the soil samples 
obtained fi-om the lateral moraines at Tsomoriri-Plains. Cell den- 
sities of methanogens (as determined from the mcrA gene copies) 
were close to the detection limit of 5 x 10^ gene copies per gram 
dry weight and rang ed from 5 x 10^ to 1.5 x lO" (Figure 5). The 
cell densities of cyanobacteria (estimated from the 16S rRNA 
gene copy numbers) were much higher and ranged from 2 x 
10'' to 4.2 X 10^. The cell densities of both cyanobacteria and 
methanogens were relatively constant with respect to the glacier 
moraine, with no apparent pattern. However, in all cases the cell 
densities of both cyanobacteria and methanogens were higher in 
biocrust samples compared to bare soil samples. 

DISCUSSION 

The occurrence of methanogens in dry soils of hot deserts has 
previously been demonstrated (Angel et al, 2011, 2012). Here 
we tested for their presence in high-altitude cold deserts, in both 
mature as well as in very young soUs, recently exposed by glacial 
recession. East Ladakh in the north of India is characterized by 
both very high altitudes (leading to a very cold climate) as well 
as little precipitation, resulting from its geographical location in 



the monsoon rain shadow. These characteristics together with its 
pristine nature made it an ideal location for performing these 
tests. 

Overall, aU three sites examined in this study (Nubra Valley 
glaciers, Tsomoriri Plains glaciers, and vegetation belts at 
Tsomoriri Plains) showed methanogenic potential, though not 
in all samples and to varying degrees. As expected, only some 
of the samples from the younger moraine soils produced CH4 
in our incubations, and the production rates were nearly always 
lower than in the samples from the various vegetation belts, 
where the soil is more mature. The production of endogenous 
CH4 under anoxic slurry conditions in soils obtained from high 
altitudes stands in contrast to a previous report by West and 
Schmidt (2002), who could stimulate methanogenesis in alpine 
soUs only by supplying H2 to their slurries. Hence, it seems that 
methanogens were present in their soils, but were limited by the 
supply of substrate. In the Himalayan samples, on the other hand, 
the entire microbial cascade involved in degradation of organic 
matter to CH4 plus CO2 was apparently present and could be 
activated upon anoxic incubation. 

The vegetation belt samples produced CH4 at rates reach- 
ing up to 1200 nmol gdw^^ d^^ while the moraine soil samples 
reached at most 850 and 142 nmol gdw^^ d^' CH4 for the lat- 
eral moraines from Tsomoriri Plains and the frontal moraines 
in Nubra Valley, respectively. These rates (particularly for the 
vegetation belt samples) are well within the range of values 



Frontiers in Microbiology | Terrestrial Microbiology 



December 2013 | Volume 4 | Article 359 | 8 



Aschenbach et al. 



Methanogens in high-altitude cold desert soils 




c*^ - 

<r> <? „o 

Cyanobacteria in biocrust 
Methanogens in biocrust 



cfi^ c?^ J" 
<^ ■<? .0 

■ Cyanobacteria in bare soil 

■ Methanogens in bare soil 



FIGURE 5 I Gene copy numbers of the 16S rRNA gene of cyanobacteria 
and the mcrA gene of methanogens in the biocrusts and bare soils of 
the lateral moraines in Tsomoriri Plains. Horizontal line: detection limit 
(5 X 10^ gene copies gdw"^ ). Means ± SE of n = 6 (2 independent qPCR 
reactions containing 3 replicates of the PGR mixture). 



previously reported by Angel et al. (2012) for upland soils from 
various parts the world under the same incubation conditions, 
indicating little effect of the particular local stress factors (e.g., low 
temperatures) on the potential methanogenic activity of the soils. 
Also the fact that of the three vegetation belts tested, the subnival 
one, which is the coldest but the wettest one, showed the highest 
methanogenic rates indicates that even in this cold-desert, water 
scarcity is more limiting than temperature. Similar to previous 
reports by Angel et al. (2012), biocrust samples always exhibited 
higher CH4 production rates than their bare soil counterparts. 

Analysis of the soil organic matter and nitrogen contents in 
the samples showed very low values and overall little differ- 
ence between sites (although they were generally lower in the 
Tsomoriri vegetation belts samples). Lack of substrate can there- 
fore not explain the difference in CH4 production rates between 
the samples and the lack of methanogenic potential in many of 
the moraine soils. Instead we suggest that it is the development 
of the microbial community (and particularly the development 
of mature biocrusts) that plays a crucial role in determining 
whether CH4 will be produced or not and to what extent, i.e., how 
much time has passed from the glacier recession. The observed 
CH4 production gradient along the lateral moraines of Tsomoriri 
Plains and the frontal moraines of Nubra Valley would sustain 
this theory because the oldest soils (TS3/N3) showed the highest 
rates and the greatest proportion of active methanogenic samples. 
Nemergut et al. (2007) who studied the development of micro- 
bial communities along a receding glacier in Peru reported not 
just an increase in the complexity of the microbial communities 
with soil age but also an overall increase in similarity between 
the communities in the young soil to that of a mature soil with 
age. Similarly, Gangwar et al. (2009) could show an increase in 
both bacterial population size and diversity with increasing alti- 
tude in Western Himalayian mountain tops. A similar trend was 
reported for the Rocky Mountains (Bryant et al, 2008). This 
indicates that despite being considered "dispersion unlimited" 
(Lindstrom and Langenheder, 2012) soil bacterial communities. 



much like plants, undergo successional development over time 
and do not simply grow in size from the moment of soil expo- 
sure. 

Analysis of the stable isotopes signatures as a reflection of 
the pathways leading to the formation of CH4 showed a clear 
dominance of the hydrogenotrophic pathway (~50-100%) dur- 
ing the first third of the incubation period, but a shift toward 
the aceticlastic pathway in the latter part of the incubation. While 
this is a typical pattern for such slurry incubations and has often 
been reported in the past, also for upland soils (e.g., Roy et al., 
1997; Glissmann and Conrad, 2002; Angel et al., 2012), this might 
not reflect well the processes as they occur in the field, partic- 
ularly for aerated soils. Angel et al. (2012) showed that in the 
presence of oxygen acetate is consumed by other processes (prob- 
ably heterotrophic ones) and CH4 is produced exclusively from 
H2/CO2. 

Initiation of CH4 production commenced after a lag of 19-34 
days on average. The lag time for the mature soils from the 
various vegetation belts was in agreement with lag times mea- 
sured for other upland soils, but the moraine soils needed sig- 
nificantly more time to initiate CH4 production (Angel et al., 
2012). Assuming no inhibitors in the soil and accounting for 
the very low concentrations of alternative electron acceptors 
such as NO^^ (>2mg kg^'; lanatkova et al., 2013) we assumed 
that the long lag times resulted from the low population size 
of the methanogens. Indeed, the initial population density of 
methanogens in the soil from Tsomoriri moraines was very 
low (lO^-lO'* copies gdw ^) and in several cases even below 
the detection limit. However, it was about tenfold larger in 
the biocrusts, showing again that biocrusts are favorable habi- 
tats for methanogens despite their photosynthetic nature. StUl, 
these values are lower than those obtained for mature biocrusts 
from the Negev Desert (Angel et al, 2012), indicating either 
that these Himalayan crusts still undergo development, or that 
cold desert crust inherently harbor a smaller methanogenic 
community. The numbers of methanogens were also much 
lower than the densities observed in anoxic soils such as rice 
fields (usually 10^ to 10'' copies gdw^'; Conrad and Frenzel, 
2002). 

Cyanobacteria are typically the dominant primary producers 
in biocrusts and also act as the pioneer colonizers in the suc- 
cessional development of the crust (Biidel, 2003). In a previous 
study it was shown that cyanobacteria comprised 70-99% of all 
phototrophs in these Himalayan slopes, and were particularly 
dominant in the newly exposed moraines (Rehakova et al., 201 1; 
Janatkova et al, 2013). Also our quantification of cyanobacteria 
in the samples from Tsomoriri Plains showed a high density of 
these organisms even in these young soils. Nevertheless, num- 
bers of cyanobacteria were lower in the bare soil samples than 
the biocrusts and were overall lower by about an order of mag- 
nitude compared to mature crusts from hot deserts (Steven 
et al., 2012; Angel and Conrad, 2013). The association between 
cyanobacteria and methanogens is well-known in marine mats 
(Hoehler et al, 2001) and has also been postulated for desert 
crusts (Angel et al., 2011). This notion is strengthened here with 
the occurrence of methanogens in these newly developed crusts 
alongside cyanobacteria, which are typical for mature crusts 
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(e.g., lichens and mosses; Belnap and Eldridge, 2003). As the 
main nitrogen fixers in arid environments (Evans and Ehleringer, 
1993) cyanobacteria are crucial for the proliferation of nearly 
all other microorganisms in the soil, including methanogens. 
Brankatschk et al. (2011) pointed out that N2 fixation, mineral- 
ization, nitrification and denitrification are important drivers of 
N turnover in young soils and that soils along a glacier forefield 
were characterized by a high abundance of N2 fixing organ- 
isms. Also Nemergut et al. (2007) showed that cyanobacteria 
played an important role in the soil development of unvege- 
tated, recently deglaciated soil and that the cyanobacteria were 
the abundant N2 -fixing clade. Apart from supplying nitrogen it is 
also possible that methanogens in biocrusts rely on cyanobacteria 
for the supply of degradable organic matter or even hydro- 
gen, as is known to occur in hypersaUne mats (Hoehler et al, 
2001). 

Fingerprinting of the archaeal community in the native 
soils showed that they were dominated by Thaumarchaeota, 
which is common for upland soils (formerly considered to be 
Crenarchaeota; Bates et al, 2011), and even for desert biocrusts 
specifically (Soule et al., 2009; Angel et al., 2012). Only few mem- 
bers of the thaumarchaeotal phylum have been cultivated so far, 
all of which seem to belong to the ammonia oxidizing archaea 
guild (Offre et al., 2013). Yet considering the abundance of 
Thaumarchaeota in natural environments it cannot be excluded 
that they perform other biogeochemical processes than ammo- 
nia oxidation. The dominance of Thaumarchaeota indicates that 
the methanogenic population comprised less than 1% of the ini- 
tial archaeal population (the typical detection limit of TRFLP). 
In fact, considering our quantitative measurements of ~10^ 
methanogens gdw^^ and an average population size of ~10^-10^ 
copies gdw^^ for archaea in desert biocrusts (Soule et al., 2009; 
Angel and Conrad, 2013) methanogens most likely comprised 
less than 0.01% of the archaeal population. This estimation of 
total archaeal population size is also supported by recent quan- 
tification of soils from this environment (RA, unpublished). 
Similar to previous reports on methanogens from upland soils 
(both arid and humid), this high-altitude cold desert was also 
dominated by methanogens of the genus Methanosarcina and 
Methanocella (Nicol et al., 2003; Poplawski et al, 2007; Angel 
et al, 2012; Scavino et al., 2013). This is in stark contrast to 
the communities found in typical methanogenic environments, 
such as wetlands, rice fields, lake sediments, and guts of rumi- 
nants or termites, which host a variety methanogenic genera 
simultaneously (Lueders et al., 2001; Sjoeling and Cowan, 2003; 
Banning et al., 2005; Brune, 2011; Nicholson et al., 2007). It has 
been postulated that the occurrence of only these two genera of 
methanogens is related to their relative oxygen tolerance com- 
pared to other methanogens (Erkel et al., 2006; Angel et al., 2011), 
and the findings of this study further corroborate the notion that 
Methanosarcina and Methanocella are globally distributed upland 
soil methanogens. However, the environment in these Himalayan 
mountaintops poses yet another type of stress in the form of 
extreme cold temperatures and frequent freeze-thaw cycles. It 
is therefore not surprising that two previous studies found a 
majority of psychrophilic and psychrotolerant bacteria in the 
communities from these soils (Nemergut et al., 2007; Gangwar 



et al., 2009). Hence, Methanosarcina and Methanocella appear to 
be also psychrotolerant in addition to their ability to tolerate air 
exposure and desiccation. 
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FIGURE Al I (A) Schematic depiction of the location of sampling sites 
along a glacier and code of sample denotation. (B) Overview of soil samples 
taken at different positions at Tsomoriri and Nubra. In Tsomoriri samples of 



matured soil were taken from three vegetation belts and of young soil from 
the lateral moraines of a glacier. In Nubra young soil samples were collected 
from the frontal moraines of three different glaciers. 
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FIGURE A2 I Archaeal 16S rRNA phylogenetic tree. Maximum likelihood 
phylogenetic tree based on aligned partial 16S rRNA gene sequences obtained 
from clones from four incubated soil samples. Clusters containing sequences 



obtained in this study are shown in colors. The length of the terminal 
restriction fragment is indicated in brackets. The number of sequences 
obtained is indicated in each cluster. Scale bar: substitutions per base. 



www.frontiersin.org 



December 2013 1 Volume 4 | Article 359 | 13 



Aschenbach et al. 



Methanogens in high-altitude cold desert soils 



Table A1 | Methane production rates and lag times of all samples. 
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"Means of the 3 technical replicates. 

^Samples from Nubra Valley were all biocrust samples. 

"BLD, below the limit of detection. 
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